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ABSTRACT 

This paper examines the possibihty that hadronic processes produce the >100 MeV 
photons in the prompt phase of gamma-ray bursts (GRBs) observed by the Fermi- 
LAT. We calculate analytically the radiation from protons and from secondary electron- 
positron pairs produced by high energy protons interacting with gamma-rays inside of 
the GRB jet. We consider both photo-pion and Bethe-Heitler pair production processes 
to create secondary electrons and positrons that then radiate via inverse Compton and 
synchrotron processes. We also consider synchrotron radiation from the protons them- 
selves. We calculate the necessary energy in protons to produce typical Fermi-h AT 
fluxes of a few /zJy at 100 MeV. For both of the photo-pion and Bethe-Heitler pro- 
cesses, we find that the required energy in protons is larger than the observed gamma- 
ray energy by a factor of a thousand or more. For proton synchrotron, the protons 
have a minimum Lorentz factor 2 x 10^. This is much larger than expected if the 
protons are accelerated by relativistic collisionless shocks in GRBs. We also provide es- 
timates of neutrino fluxes expected from photo-hadronic processes. Although the flux 
from a single burst is below IceCube detection limits, it may be possible to rule out 
photo-hadronic models by adding up the contribution of several bursts. Therefore, 
photo-hadronic processes seem an unlikely candidate for producing the Fermi-hAT 
radiation during the prompt phase of GRBs. 

Key words: gamma rays: Bursts - radiation mechanisms: non-thermal 



1 INTRODUCTION 

The Fermi satellit e has detected 29 g amma-ray bursts (GRBs) with photons of energies >100 MeV with the Large Area 
Telescope (LAT) l Atwood et al.ll2009l ). The photons detected during the prompt phase of GRBs by the Fermi- LAT are 
emitted with a delay of a few seconds for long GRBs, and are observed for a longer du ration of time comp ared to the lower- 
energy photons (~1 MeV) observed by the Fermi Gamma-ray Burst Monitor (GBM) (jAbdo et al.ll2009br ). Hadronic models 
could explain the delayed emission detected by the Fermi-hAT because of the additional time needed to energize the protons 
( Razzaaue et al. 2O10l : Asano fc Meszar"o^ 2012 1. Since protons do not suffer from radiative losses like electrons, it is possible 



that GRBs accelerate pro tons much more efficiently, causing the hadronic emission to dominate at high energies, (jAsano et al 



20091 : iMurase et al]|2012l ). The relative efficiency of accelerating protons or electrons in GRBs depends on the mechanism 
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Blandford fc Eichlerlll987r ) or magnetic dissipation {e.g ^ . . _ _ . , _ . .. _ _ . 

Lvutikovll2006 ). In this work, instead of dealing with the detailed specifics of particle acceleration 
for realistic GRB jet models, we assume that high energy protons exist as a power-law distribution and calculate the energy 
required in protons to reproduce the high-energy photon flux observed in the LAT band during the prompt phase. 

There is a growing body of evidence that the observed emission in the LAT band after the prompt phase may be due to 
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synchrotron radiation from the ear ly afterglow emission, where electrons from the medium surrounding the GRB are acceler- 
ated by the external forward shock ( Kumar &: Barniol Durarj2009l : Gao et al]|2009l : Corsi et al ] |2O10l :l iKumar fc Barniol DuranI 
The most convincing evidence of external forward shock producing the > 100 MeV photons is that the late tim e optical 



and X-ray afterglow data can accurately predict the earlier flux in the LAT band I Kumar fc Barniol DuranI 20091 ). Under 
the external forward shock origin for the LAT emission, the delay in the LAT is due to the delayed onset of the external 
forward shock emission. The external-forward shock emission begins when the jet emitted by the central engine deposits half 
of its kinetic energy to the surrounding medium. Therefore, there remains the possibility that while the later time (T > TgcJ^) 
LAT-band flux is from the early afterglow, the earliest high energy photons observed during the prompt phase are created by a 
different process. The external forward shock is unable to explain the claimed variability present in the LAT prompt emission 
i Maxham et al.ll201ll ). Furthermore, although the most G RB spectra can b e fit with a smoothly-joined broken power law 
extending several decades in energy (i.e. the Band function, Band et al. 19931 ). a few bursts exhibit deviations from the simple 
Band function. Notably , GRBs 090902B, 090926A, and 090510 sh ow evidence of an additional high-energy spectral compo- 
nent (jZhang et al.ll201ll : lAbdo et al.ll2009al : lAckermann et al.ll2011^ ■ This spectral feature is usually transient and disappears 
before the prompt phase is over. For example, in GRB 090510, which has a Tqo of 1.5 s, it be comes statistically impossible to 
distinguish the additional power-law component ~ 0.8 s after the trigger (|Asano et al.ll2009l ). Similar spec tral evolution was 
seen in GRB 090926A, where an additional power law component appears after ~ 11 s, and Tgo is 13.6 s (jAckermann et al. 



201ll ). This suggests that the high-energy prompt emission may have a different origin than the extended high-energy emission 



for some GRBs. 

This paper is similar to Asano et al. ( 2009[ ) and Asano et al. ( 2O10l l in that it addresses interactions between high-energy 
protons and the low-energy (GBM) photons as a possible mechanism to produced the high energy photons observed in 
_Fermj-LAT detected GRBs, but it differs in a few major ways. First, this paper investigates hadronic models of radiation — 
i.e. photo-pion production, Bethe-Heitler pair production, and proton synchrotron not only to explain the extra spectral 
components of the few GRB that exhibit them, but also to investigate the possibility that hadronic models are responsible 
for the majority of the GRB prompt emission observed in the Fermi-hKT. Second, instead of using Monte Carlo methods or 
detailed numerical codes, we use simplifying assumptions to calculate everything analytically or semi-analytically whenever 
feasible. Finally, the most important point of this paper is to provide analytical estimates of the energy required in protons to 
explain the LAT flux with hadronic emission. These estimates allow us to see how the energy requirement depends on various 
parameters. 

The LAT and GBM detected photons are likely to have separate origins during the prompt phase because of the measured 
delay between them. Since we are interested in explaining the observed LAT flux of GRBs, we will assume that the < 100 
MeV portion of the GRB emission is of unknown origin and perfectly described by a Band function flt with an observed peak 
energy of i/p ~ 1 MeV (the exact value, is of course, different for different bursts, and Vp also evolves with time during a single 
burst). The Band function has a low energy index a and a high energy index /?. We only consider interactions between the 
photons from the Band function and high-energy protons in the GRB jet. That is, we do not consider second-order processes, 
such as proton colliding with a pion-produced positron, et cetera. This choice greatly decreases the computational difficulty 
but produces results that only hold in regions where the contributions from these processes are likely small, i.e. above the 
photosphere. In two zone models, the l ower energy emis si on is produced at a si gnificantly smaller radius than the radius where 
high energy emission is produced. As Zou et al. ( 201lh : Hascoet et al. ( 20121 ) have noted, a two zone mode l will reduce the 
minimum Lorentz factor by a factor of ~ 2 — 5 com pared to the value obtained using a one zone model such as lLithwick fc Sari 
l 200ll ): lAb"do et all |2009bh : iGreiner et all (|2009h . A two zone model is unlikely to change our results by a significant factor. 
We provide the dependencies of the efficiency on Lorentz factor to help determine if a lower Lorentz factor could make the 
energy requirement more attractive. 

The paper is organized as follows: In section [2] we calculate the positrons produced through the photo-pion process. In 
section 12.11 we provide an analytical estimate for the required luminosity in protons to match a flux of 1 /xjy at 100 MeV 
with the photo-pion process and provide the corresponding neutrino flux in l2.1.1l The required proton luminosity calculation 
is then repeated in section 12.21 however it is done numerically and more accurately. In addition, the maximum efficiency of 
the photo-pion process is calculated for given GRB parameters. In section [3] we calculate the electrons produced through 
Bethe-Heitler pair production and compare this to the photo-pion calculation in W2.2\ Section |4] contains a calculation of the 
maximum efficiency for the proton synchrotron production of typical fluxes for LAT detected GRBs, as well as an expected 
neutrino flux from the proton synchrotron model. Finally, in section [5] we summarize and discuss our results. For calculations 
of the luminosity distance we take cosmological parameters of Ho — 71 km/s/Mpc, Sim. = 0.27, SIa = 0.73. 



^ Tgo is the time at which the GRB has radiated 90% of its observed photon energy 
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2 PHOTO-PION PRODUCTION 

The photo-pion process refers to the production of pions (tt'' and tt*) in collisions of photons with protons. The decay of these 
pions produces high energy electrons and positrons that can then produce high energy photons via synchrotron radiation. 
The photo-pion process is likely to be important in situations where electrons are unable to be accelerated efficiently to very 
high Lorentz factors, but protons are. It also offers a way to beat the well known limit on the maximum synchrotron photon 
energy of ^ ^O j^ M eV for shock-accelerated electrons, where F is the Lorentz factor of the source and z is the redshift 



I de Jager et al. 



19961 ). 



The delta resonance, p"^ 4- 7 — >■ , has the largest cross section and has the lowest energy threshold of the photo-proton 
resonances and is therefore the most important photo-pion interaction to consider. The delta resonance has two main decay 
channels. A"*" — >■ tt"^ + n and A^ — 7> tt" The neutral pions quickly decay further as tt" — >■ 7 + 7. The threshold photon 

energy for photo-pion production is approximately 200 MeV in the proton rest frame. The photon energy at the peak of the 
GRB spectrum, in jet comoving frame, is fpll + z)/T, where the peak frequency in the observer frame is Up and the GRB jet 
is moving with a Lorentz factor P. For a proton to undergo pion production with photons of energy Vp, the proton Lorentz 
factor must satisfy 

7p>2x 10*F2/.p-^(l-K2)-\ (1) 

where Vp^ is the observed peak frequency in units of 10® eV; the standard notation X„ = is used. At the threshold, the tt" 
is produced more or less at rest with respect to the proton. This means that the tt" will decay into two photons with energies 
~ 100 F2i'~g(l -I- z)~^ TeV in the observer rest frame. This energy is well outside the Fermi-'LNJ: band. The high-energy 
photons created through 7r° decay will interact with the lower energy photons to produce pairs. If the optical depth of the 
7 + 7 pair production for 1 TeV photons is much greater than one, they will readily produce pairs with Lorentz factors of 
~ 10" F2t'~g(l -|- z)~^ , a similar value to the electrons produced by tt^ decay, see eq ([2]). If the optical depth of 7 -I- 7 pair 
production is much less than one, then the photons will escape the GRB jet, and the vr'' will not affect the ~ 100 MeV flux. 

The tt"*" decays as 7r+ — >■ /i+ -\-Vp,, and the anti-muon decays further as fj,^ — >■ e+ + + Isospin conservation arguments 
suggest a branching ratio for the delta resonance for tt"*" : tt^ of 1 : 2. However, for protons interacting with a power-law 
distribution of photons, where there are a sufficient number of high-energy photons to excite higher energy resonances as well 
as allow direct pion production, the ratio of charged pions, tt^, to n eutral pions is actually closer to 2 : 1. This ratio is more 



or less independent of the photon index ( Rachen fc Meszaros 19981 ). As an approximation, we take the cross section of the 



delta resonance, but only consider the high energy electrons produced by the tt* decay. This underestimates the high energy 
electron production rate by a factor of 3 when the GRB jet is opaque to photons from the vr" decay. 



2.1 Analytical Estimate 

As with the vr", at the threshold energy, the tt^ (and subsequently the are produced almost at rest in the rest frame 
of the proton and therefore have the same Lorentz factor as the proton given in eq ([T]). On average, the positron carries 
roughly one-third of the energy of the muon (the remaining two-thirds goes to neutrinos). The Lorentz factor of the in the 
comoving jet rest frame is 

7e~i^7P~10'r2^,^^(l + ^)-\ (2) 
and me are the muon and electron masses respectively. 

By requiring that a typical positron produced through pion decay, with Lorentz factor given by eq ([2| and charge q, 
radiates at a desired frequency v (~ 100 MeV for Fermi-LAT), we solve for the magnetic field in the jet rest frame, B, in 
Gauss. 

IEIeL — ^ ^ 1.6 X IQ-^i/g erg ^ B ~ imusul Jl + zfv:;^ Gauss (3) 

27rmeC(l + z) 

This value for B requires the minimal proton energy to match an observed flux. If the energy requirements are too large when 
B is equal to eq (O, they will be even worse when the magnetic field is not equal to eq (12|g 
The observed specific synchrotron flux, is 

^/?,q''BN,T{l + z) 



■^^ = 4^dim> ^ iV,,5oBFd-48(l + (4) 



Equation ((Sjl is close to the magnetic field value requiring the minimal proton energy when the peak frequency of a positron with 
Lorentz factor given by eq l(2| is above the cooling frequency i^c. If a typical photo-pion produced positron is not cooled by synchrotron 
radiation, for typical GRB spectra and efficient proton acceleration, the minimum necessary proton luminosity will occur at a B such 
that Uc is 100 MeV. For the vast majority of allowed GRB parameter space, the photo-pion peak is cooled, so the statement that eq ((Sjl 
is a best case scenario holds. 
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where A^e is the number of electrons that radiate at v, and cLl is the luminosity distance to the GRB. Thus, the number of 
e"*" needed to produce the observed flux at v is 

iV..8xlO-|g^. (5) 

where /i^.^jy is the observed flux in /iJy. 

The number of protons, with energy above the pion production threshold, required to produce the necessary electrons in 
eq ([5]) is Np ~ ^Ne/rp^, where Tp^ is the probability that a photon of frequency ~ !/'p(l + z)/r collides with a proton with 
a Lorentz factor given by eq ((T)) and produces a pion. This probability is approximately equal to the optical depth to pion 
production and is given by R/r, where apj is the cross section for the delta resonance, ap-, = 5 x 10 cm . 

n-y is the number density of photons in the comoving frame of the GRB jet, and R is the distance from the center of the 
explosion. For a GRB of observed isotropic luminosity (integrated over the sub-MeV part of the Band function spectrum) 
and observed peak frequency Vp (in eV), the number density of photons in the comoving frame of a GRB jet is 

"- = 4vrCTc.,(1.6xlO-i^erg/eV) "'"^° ^..52i^i5 .pA^ + z) cm. (6) 
This gives an optical depth of 

Tp^ ^ 0.8L^,52RT5^2'''^pM^ + ^y'- (7) 

Using ((6| and (O the number of protons needed to produce a specific flux, f^, is 

Np ^ 10'**/^,^JyCii,28r2-Rl5!^p,6-B ""^I/^ 52- (8) 

The corresponding energy in these protons is 

r, r on , inSl r^/i/.nJy'^i, 28-^15 

Ep ^ Npi-YpiripC )r ^ 3.0 x 10 ' erg. (9) 

-D-tj-y, 52(1 + Z) 

It is more useful to consider the luminosity carried by these protons, Lp, as this can be directly compared to the observed 
7-ray luminosity, L^. The ratio L^/ Lp will allow us to determine the efficiency of the photo-pion process for the generation 
of > 100 MeV 7-rays. The proton luminosity is related to Ep by 

Lp = EpT X max { t^y^ , t^^^^ } , (10) 
where fdyn is the dynamical time in the jet comoving frame, 

tdyn = nOR^sV^^ S, (11) 

and tcooi ~ JemeC^ / Psyn is the synchrotron cooling time in the jet frame. Psyn is the total synchrotron power radiated by a 
positron. Using the magnetic field in the jet comoving frame given in eq ([31 , tcooi becomes 

= « 8 X 10- 2 ^5 - (12) 

Inverse Compton cooling is neglected because the electrons considered here have a Lorentz factor of > 10®; so the IC radiation 
is greatly decreased because of Klein-Nishina suppression. 

Substituting Equations ((3]), ((9l), pip . & (|12p into (|10p . we find the minimum proton luminosity necessary to match a flux 
of fv,njy at vs with the photo-pion process is 

'2 X 10^^ riL-^52'^8"'^"6/-,MJyrfi,28(l + erg S'^ fdyn < tcool 

Lp>{ ' ' (13) 

4 X 10^^ ri_Ri5L~52i/8t'p,6/^,MJydi,28(l + erg S"^ t^ool < tdyn- 

At first glance, the proton luminosity does not seem prohibitively large. However, the strong dependence on F has the potential 
to increase the proton energy requirement tremendously. 

To assess the viability of the photo-pio n process producing the > 100 MeV photons detected by Fermi, let us consider 
a bright Fermi-LAT burst, GRB 080916C jAbdo et al.lboOQbl '). This burst was detected at a redshift of 4.3, which has a 
corresponding dL_28 ~ 12. The peak of the observed spectrum was at 400 keV, and the fiux at 100 MeV during the prompt 
emission was /i, ~ 3 fJ.3y- The 7-ra y isotropic l uminosity for GRB 080916C was 1/7,52 ~ 20, and the minimum jet Lorentz 
factor was estimated to be F2 ~ 9 ( Abdo et al.ll2009b \ For these parameters, we find tcooi < tdyn as long as i? > 10^^ cm, 
implying that the required luminosity in protons with 7^ > 10^ is Lp --^ 1.5 x lO'''" -R15 erg s~^. This is a factor of ~ 700 
times larger than the 7-ray luminosity at i? = 10^^ cm. Below R — 10^'', tdyn < tcooi and the proton luminosity has no R 
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depe ndencelfl This efBciency is below the efficiencies of the order 10-2 0% estimated for other GRBs using afterglow modeling 
fe.o JPanaitescu fc Kumail2002l : Fan fc Piranlbood : IZhang et al ] |2007h and makes the photo-pion process an unlikely candidate 



to produce the prompt-LAT emission observed in this GRB. 



2.1.1 Neutrino Flux 



In addition to producing high energy electrons, photo-pion production also results in high energy neutrinos. For the photo-pion 
production models of the observed LAT emission, it is possible to directly correlate the flux at 100 MeV to an expected flux 
of neutrinos. Since there are two muon neutrinos created for every in tt^ decay, we can simply use eq ^ and ([3} to find 
corresponding number of neutrinos, 

= m"^ .f..f.iydl2sUs'u~lrUi + (14) 

These neutrinos will have the same energy as the electrons on average, with an observed energy of 

Ey = lifvlv-lil + zY^ GeV (15) 
This corresponds to an observed flux, Fy, at 10^ + ^Y"^ GeV of 

Fy = Ey{l + z)NyV X ^ max Ci}, (16) 

( 5 X 10-^ /^,MJy^r5'r|l'8"V-^(l + ZY^ GeV Cm'^s'^ fdyn < tcool 

Fy^l ' (17) 

[lO ^ fv.fiJyfH GeV cm ^S ^ tcool < tdyn. 

The neutrino flux does not depend on the luminosity of the GRB, as it is fixed by the flux at i/g, the peak of the photo-pion 
synchrotron emission. When the photo-pion electrons are in the cooled regime — as is expected for much of the GRB parameter 
space — the neutrino flux does not depend on any of the jet parameters. The neutrino flux depends only on the observed LAT 
flux at 100 MeV; when the electrons are in the fast cooling regime, the synchrotron flux depends directly on the electron 
energy flux, and the neutrino flux depends directly on the electron energy flux. Of course, the neutrino energy does depend 
on r, as seen from eq (|15p . 

To get a rough idea of an expected neutrino count-rate at IceCube, we fit the averaged effective area for muon neutrinos 
at IceCube given in lAbbasi et al] ||2012| ') by A ^ 100 m^ x ^Ey/{IQQ TeV). Our fit agrees well with the averaged effective 



area of 59-string detector when E^ > Sx 10* GeV. We expect the following neutrino counts per second of LAT emission from 
the photo-pion process: 

/i^.^Jy-Rrs^rli/g'Vp g'^(l + z)^* counts S"^ tdyn < icool 

(18) 

10"^ fu,t,jyT2 Vgi^p.^l + z) counts s"^ t 

cool ^ ^dyn- 



dt 



Therefore for a bright GRB detected by the Fermi-LAT with fv,ij.jy ~ 2, r2 ~ 9, 2: ~ 2, Vp^, ~ 1, and T90 = 10 s. We find 
that tcooi < tdyn if -R > 10^^ Cm. If i? > 10^* cm, we expect about 0.07 neutrinos of energy 9.0 x 10^ GeV over the course of 
the burst. If 7? < 10^^ cm, the number of neutrinos is increased by a factor R^r^ and the energy stays 9.0 x 10^ GeV. And so 
if we add up the contributions from all Ferrrai-LAT bursts, we expect IceCube to detect ~ 1 neutrino. 

2.2 Numerical Calculation 

The calculation presented in i]2.1l is a rough estimate to the energy requirement for protons, but it doesn't give any spectral 
information about the photo-pion radiation and assumes that all the protons have the same energy. Furthermore, it doesn't 
take the finite width of the delta resonance into account. All of these corrections go in the direction of decreasing the efficiency 
of the photo-pion process. A more rigorous calculation is presented in this subsection. 

The distribution of photons in the jet rest frame is assumed to be isotropic, and the Band function is approximated by 

J(^) = n., x<( (19) 



^ See previous footnote about how the luminosity estimate may be too pessimistic when tdyn < tcool- A more accurate calculation that 
accurately minimizes Lp with respect to B for a given R, F, etc, is presented in the i|2.2l 




6 P. Crumley and P. Kumar 



Here ep — hvp{l + z)/{rmec'^) ~ .02r^^i'p,6(l + z) is the dimensionless photon energy at the peak of the spectrum in the 
jet comoving frame, and emin is the dimensionless photon energy below which the Band function no longer fits the observed 
GRB spectrum, emin is poorly constrained by GRB observations, and depends on the prompt radiation mechanism and radius 
of emission for the prompt emission. If emin corresponds to the synchrotron self-absorption frequency, it is likely of order 
10"^, corresponding to an observed frequency of a few eV. We more conservatively assign an emin value based on the lowest 
observed frequencies by Fermi Umin ~ 1 keV, with a corresponding emin ~ 2 x 10~^r^^j/min,3(l + z). In reality, the choice of 
emin has very little effect on the synchrotron flux in the Fermi-LAT band, as the photo-pion positrons and electrons produced 
from a proton interacting with photons of energy emin will be of very high Lorentz factor, 7e ~ 10^. For our calculation, emin 
corresponds to a v^iin of ~ 1 keV and is included in our calculations only for completeness. 

If Q ~ 1 and the GRB has an isotropic luminosity Lj, then the value for the number density of photons per rrieC^ at ep. 



Wej, = 8 X 10^^L^,52i?i5^fc'p,6 ^ (1 + z) ^ cm ^. 

We assume that the proton number distribution is a power law. 



djp Ji <Jp < 7n 



(20) 



(21) 



with a minimum Lorentz factor ji ~ 10 and a maximum Lorent z factor give n by requiring the protons to be conflned to the 
jet, i.e. the Hillas criterion, 7max = rm^i = 3 x lO^'-B-RisF^^ ( HillasI 1984 ). Np^i is the number of protons in the emitting 
region of a GRB between 7^ and 7; + d7p. As these high energy protons travel through the jet, they will interact with the 
photons that make up the Band function, creating secondary particles. The total interaction rate, Np.y, for a proton with 
Lorentz factor 7p and a photon with energy e depends on the angle- integrated cross section: crp-y(e'). e' is the energy of the 
photon in the nuclear rest frame, e' = 7pe(l — Ppn), where is the cosine of the angle between the proton and photon and /3p 
is the velocity of the proton divided by c. The interaction rate is 



Np. 



c 

4-K 



dQ, / de n{e,Q,){l — I3pfi)apj{e'). 



(22) 



Since 7p ^ 1, /?p ~ 1, and we are approximating the sub-MeV Band photons as isotropic in the rest frame of the jet, the 
number of scatterings is approximated by 



dNp 



c dN„ 



de 



de e'(7p^(e'). 



(23) 



d7p 27^ d7p . 

We approximate the cross section of the delta resonance, (7p^(e'), as 5 x 10~^* cm^ if 530 < e' < 760 and otherwise. As before 



we treat the pion and muon as decaying instantaneously without any energy losses and approximate 7e 
l|23[) . re- written in terms of the produced electrons, is 



d'ye 



c dNp djp 



, n{e) , , , , 

de — — I de e apj[e ) 

^ Jo 



707p. Equation 



(24) 



lO^ji d-yp d'ye jg 

When evaluating the scattering rate, it is convenient to define two electron Lorentz factors of interest: 7peak, the electron that 
is produced from a proton interacting with photons at the observed peak in the gamma-rays, and 7brGak, the electron that is 
produced from a proton interacting with the lowest energy photon in the Band function, which is taken to be fmin- 7pGak and 
7brcak are equal to 

7peak = 1.2 X 10''r2t'p^6(l + 2)"' (25) 
7tacak = 1.2 X 10'r2//min,3(l + 2)"'- (26) 

Carrying out the integration in equation (|24p . we find the rate of electrons produced through the photo-pion processes is 

/ 7, y-"-' 

Ne.r, I ) 707i ^ 7e < 7pcak 



dN 
d'ye 



Ne. 



Ne. 



0—p—l / \ Q—p — 1 

'Ypcak \ I '^fe 



707, 



'Tpcak 



8~p—l/ \ a — p—1 / \ —p — 2 

7peak \ / 7brcak \ / 7e 



707^ 



7pcak 



7break 



7pcak < 7e ^ 7broak 



7brcak < fe 



Ne. 



0.3(6 X 10-^)''A^p,,7f-'i7,52i?l5'<6' (1 + zf-^r^^/(p + 1). 



(27) 



(28) 
(29) 
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We now derive the previous estimate of how much energy the protons would need to carry to produce the observed Fermi-LAT 
flux at 100 MeV. 

For simplicity, we set a and /3 to the typical GRB parameters a — 1 and j3 = 2.2. We assume the protons have a power 
law index of p = 2 and ji — 10, corresponding to efficient acceleration in shocks. Then, the majority of the energy in the 
photo-pion electrons is contained in the electrons with 7pcak ^ 7e ^ 7brcak- This section of the power law is 

^ = 7 X W~''N,,,L,,s2RT,''^-lil + zy'T^' (i^) ' • (30) 

To calculate the total number of electrons produced, we solve the following continuity equation 
dNj-y,) d dN 

__ + _{^,iV(7e)} = ^. (31) 

We approximately solve this equation by following the standard procedure of breaking up the continuity equation into 
two regimes: one where the electrons are cooling slowly, i.e. t^yn < tcooi, and another where cooling losses are important, 
idyn > icooi- The solution of differential equation (|31|) is then approximately: 

dN 

yn ~J ^dyii ^cool 

NM = { (32) 
1 r ^ dN 

-T^ J ^dyu > icool, 

Two different cooling mechanisms are considered: synchrotron and inverse Compton cooling. To see which is the more dominant 
cooling process, we compare the power radiated by each process. The synchrotron power for electrons with 7e = 7poak is 

Psyn(7pcak) = 1.6 X 10-^ri(l + z)-'!/-^^' Crg/s. (33) 

From the condition tcooi < idyn, the electrons at the peak of the distribution will be cooled via synchrotron if 

B>2 Ri^'-''u^^i{l + z)°-'' Gauss. (34) 

This is a low value of the magnetic field, so synchrotron cooling losses are important to consider. However, for completeness, 
our estimate will consider both possibilities, when 7pcak is above cooling and below cooling. 

For inverse Compton losses, while the energy density in the photons can be very large, particularly at distances less 
than ~ 10^'' cm, the inverse Compton radiated power is greatly reduced due to Klein-Nishina suppression. For the electron 
Lorentz factor given in eq (I25|l . all of the prompt sub- MeV emission will be in the Klein-Nishin a regime if 7Deak£m in > 1, or 
i/p,6 < 201^111111,3 • The power radiated due to IC scattering in the Klein-Nishina regime is given in iBlumenthall l|l971f) : 

PKN{l) = niJnrl rfe i ^ (log (47^) - ^) . (35) 

where ro is the classical electron radius. For 7 = 7pcak, neglecting the logarithmic dependencies of variables and assuming 
Q ~ 1, the IC radiated power is 

PKiv(7pcak) = 2 X lQ-^L-,,52R^iu-l{l + z)-^ \ . (36) 



From equations p3l) and p6|) . the ratio of synchrotron power to Inverse Compton power at 7pcak is 



If we define es as the ratio of energy density in the magnetic field to energy density in radiation, the ratio becomes 

^ = 50eB,_2 { ^^^^ . (38) 

PkN \ J 

Unless es is small, the synchrotron emission will dominate over the inverse Compton emission. The inverse Compton scattered 
photons will have on average an energy in the jet's rest frame of 'ye'm^c? ~ 0.5 TeV for electrons with 7e = 7peak. These 
photons will quickly pair produce and form a cascade of secondary particles. A full treatment of this is beyond the scope of 
this paper. In any case, as can be seen in equation (|38|) . the energy in this cascade will be less than the synchrotron energy 
radiated by the photo-pion produced electrons; this allows us to ignore these pairs when for estimating the flux at 100 
MeV. 

The electron number distribution created by the photo-pion process for 7poak ^ 7e ^ 7brcak is 
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dN 
d'ye 



(39) 

5 X 10-^ B-2iVj,_,L^,52i?-2i.-l(l + Z)-Ir^^7j tdyn > tcool. 

The observed synchrotron flux, fv, at ~ 100 MeV, is calculated using the following approximation for synchrotron radiation: 

2/3 



x/3g^rBiV(7e 



ind'im.ec'^ \ 7, 



/. = (1+.)/ d-je '"i:2Z' ir {^] > (40) 



2 ^ 2nmecil + z)u 

gBr(4.13 X 10-15 eV s)' ^ ' 

Using equation ([3]) for the magnetic field value, we calculate the necessary luminosity in protons to produce 7-rays through 
the photo-pion process. The result we found for Lp is 



'3 X 10^1 riL-l2i/8-li.p-2/.,MJy4,28(l + Z)'^ erg S-^ tdyn < 

7 X 10^'' rii?i5L"52i'8i'p,6^,M.iydi,28(l + z) erg s" 



(42) 



7,52" 

In comparison to the previous estimate for Lp given in equation (|13p . the values for Lp in equation (142 |l are considerably 
larger, by a factor of ~ 100. A factor of ~ 20 is attributable to the fact that unlike the estimate given in ^2.11 this calculation 
considers protons that are part of a power law distribution that extends over several decades of energy. Additional factors come 
from the finite width of the delta resonance and keeping track of the factors that come from integration. For the parameters 
of GRB 080916C, for the expected case tcooi < tdyn, the required proton luminosity is Lp ~ 3 x 10^* iiis erg . So, the 
luminosity in the protons is 10^ times larger than luminosity in the 7-rays at R15, which is too large to be realistic for a stellar 
mass object. 

We define the efficiency, rj, as 

I 3 *-£'^,52l'8l'p,6/J^p,IyrfL,^28(l + z)"' Wn < *cool 

V^f^ = { (43) 



10 ^ r2 ^-Ri5^I/^,52l'g g/^j;jyd^_28(l + ^) ^ idyn > ^cool- 



In the previous equation, the cooled and uncooled estimates for i] were calculated by choosing a magnetic field to ensure 
the energy peak of the photo-pion-produced electrons radiated at 100 MeV. While this is convenient and pretty accurate 
maximum efficiency for analytical estimation, we also numerically calculated the maximum efficiency, allowing _B to be a free 
parameter while fixing all the other parameters {R, F, L^, Up, etc). As bounds on B, we set the minimum magnetic field 
value by requiring that the power radiated through inverse Compton is no more than 100 times the synchrotron power for 
an electron that has a synchrotron peak at 100 MeV. We set a maximum value for B such that the energy in the magnetic 
field is at most 10 times the energy in the photons. For the parameter space we considered, the B that maximized rj was well 
within these bounds. For a given L, F, R, and p, we calculate the maximum efficiency of photo-pion electrons radiating the 
desired fiux of 1 /xjy at 100 MeV. This maximum efficiency is plotted in figure [T] The part of equation (|43p corresponding to 
fast electron cooling gives an accurate prediction of the maximum ry. In the slow cooling regime, equation (|43|l predicts too 
small a value of ri; in this case the maximum efficiency is found when i? is a value such that 100 MeV is Vc- 

As illustrated in figure [T] and equation (|43|l . R, F and Lj are the only parameters capable of changing rj significantly; p 
can as well, but it is fixed by the desired photo-pion spectrum and therefore not a free parameter. From typical GRB spectra, 
we expect p to be in the rage 2.4-2.8 to match typical LAT spectra. In the bottom right panel of figure [1] we can see that p 
has almost no effect on rj when we only consider the protons creating the > 100 MeV photons (see upper red line in figure 
[1]). In figure[2l to explore how the efficiency changes with R, F, and L-^ we plotted rj in the R — V plane for various L-y. It is 
interesting to note that although rj scales as L^ for a fixed R, F, the maximum efficiency in the R — V plane only scales as 
~ L~f because the available parameter space decreases with increasing L^ due to 7 -f 7 pair production. 



3 BETHE-HEITLER PAIR PRODUCTION 

Through Bethe-Heitler pair production, protons and photons interact to create electron-position pairs directly, p + y 
p + e'^ + e~ . The Bethe-Heitler cross-section and the energy of the produced electron-positron pair depend strongly on the 
angle between the outgoing electron/positron and the proton. Therefore, it is not possible to use the integrated cross section 
to calculate the secondary electron production. Assuming that the protons and photons are isotropic in the jet's rest frame 
and using the head on approximation, i.e. the angle between the photon and proton is zero, i.e. e' — 2"fpe, the equation for 
the rate of production of secondary electrons is: 

^(7e) = 2cj^ de n{e) d^p Npi^fp)^^^. (44) 
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r p 



Figure 1. A plot of the maximum efficiencies for the photo-pion process radiating a typical _Fermj-LAT flux of 1 /xjy at 100 MeV as a 
function of R, L, F, and p. The lower blue lines corresponds to a more physically realistic case of a proton power law extending from 
a Lorentz factor of 10 to the Hillas criterion for the confinement of protons. The upper red lines are the efficiency only considering the 
energy in the protons that produce the pions and then electrons that radiate at LAT frequencies. These red lines represent an absolute 
maximum possible efficiency and corresponds roughly to our calculation in i|2.1l The dotted line corresponds to the cases when LAT 
emission could not be seen by an observer, either because the emission happens below the photosphere or because the jet would be 
opaque to photons of 10 GeV due to 7 + 7 pair production. In a two zone model for gamma-ray generation, the area where this occurs 
may differ by a small factor. When the parameters are not on the z-axis, values of = 10^^ erg/s, R = 10^^ cm, F = 800, Up = 1 MeV, 
z = 2, = 4.9 X 10^* cm, and p = 2 are taken. The photon power law indices of the Band function were set to a = 1 and /3 = 2.2. 
Since Up, a, and /3 are unable to change the maximum r] by more than an order of magnitude, their corresponding plots are omitted. 



In this equation Np{'yp) is the number of protons with Lorentz factor 7^ and n(e) is number density of photons with energy 
e. The formula for the differential Bethe- Heitler cross section , aBg , in t he Born approx imation, integrated over angles in the 
highly relativistic regime, was derived by Bethe fc Maximon (1954) (see Rachen (199^) for a more recent review). 



d7i 



3q/ctt / / 



7+ +7- + gf+T- 



lof 



27;7- 



(45) 



In this equation, 7^(7!.) is the Lorentz factor of the positron (electron), Uf is the fine structure constant, and all of the above 
quantities are in the proton rest frame. Much of the contribution to the angle-integrated cross section comes from angles 
between the photon and outgoing of order 6'^ = When 7p 3> 7±, the Lorentz factor of in the jet's rest frame is 



7± = 7p7± (1 - /3p/3± 



/ -2 

-7± 



7p 
7i' 



(46) 



Therefore, most pairs produced via the Bethe-Heitler process have Lorentz factors (in the jet comoving frame) that are smaller 
than the proton that created it. 

If e' <C nip/ me ~ 10^, the nuclear recoil of the proton can be neglected and the following equality holds: 



7-^+7- 



(47) 
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Figure 2. The maximum rj is plotted in the R — T plane for various values of to match a flux of /j, = 1 ^Jy at 100 MeV. As in 
Figure [T] we fix = 1 MeV, z = 2, dj^ = 4.9 X lO'^^ cm, and p = 2. The photon power law indices were set to typical values for the 
Band function, a = 1 and /3 = 2.2. Where 7] is not shown, the place of emission is either below the photosphere or opaque to radiation 
of 10 GeV due to 7 + 7 pair production. 



For large e , the differential cross section decreases very rapidly when 7^- < 2. Therefore, we only consider 7^- ^ 2, where the 
differential cross section is more or less constant. In this regime, the differential cross section simplifies to 

daBH 



Re- writing eq (|48 
dasH a(JT 



if 2 s; 7+ < e' 



in the jet comoving frame and using the e « 27^6, we find 



(48) 



(49) 



d7+ 2672 ' 

The integral in equation (|44|) is simplified by considering this approximate expression for the cross-section. The integral 
is now straight forward to calculate for a Band spectrum with indices a, j3 and a proton index of p. The result is 



d'ye 



/3(p+l)7. ^P^^P'* 



(7e 



lea y cr-j- Cp 
5/3(p+l) ' 



\ep-fi J V 5 



for ^ ^ 7c ^ 5/ep 



for 5/ep ^ 7e ^ 5/e 



(50) 



We now compare Bethe-Heitler pair production to the photo-pion process. The integrated cross section for Bethe-Heitler 
process is roughly 10 times larger than the cross section for the photo-pion A-resonance. For any given proton Lorentz factor 
7p, the photon energy required for Bethe-Heitler is roughly 50 times smaller than for the A-resonance. For a given 7^, the 
photo-pair will have an average Lorentz factor of ~ 7p/5 while the delta resonance will decay to a electron with an energy 
~ 707p . Consider the case where protons with a power-law distribution function with index p are scattering with a isotropic 
photon power-law spectrum n(e) cx e~". The ratio of the number of above a fixed Lorentz factor generated by Bethe-Heitler 
process compared to those generated by photo-pion process is ~ 2 x 10 x (10^)"~i x {300)~^'^^ — the first factor comes from the 
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fact that Bethe-Heitler produces a electron-positron pair compared to a single positron produced in the delta-resonance, the 
second factor is the ratio of the total cross sections for the Bethe-Heitler and photo-pion scatterings, the third factor accounts 
for the larger number of photons that participate in the Bethe-Heitler process (the threshold energy for Bethe-Heitler is 
f ~ 7(7^ and the threshold energy for photo-pion is ~ and the final factor is due to the fewer number of protons 

that can create electrons with Lorentz factor > 7^. This means that which process dominates depends strongly on which 
part of the Band function the protons are interacting with to produce with Lorentz factor > 7e. For 7e > 10^, the energy 
threshold for both processes lies below the peak of the Band function, so a = a « 1. Thus, in this regime, the photo-pion 
pairs dominate. 

However, for 7e < 10^, the threshold photon energy for both process is above the peak of the gamma ray spectrum, 
so a = /3 « 2.2 and the Bethe-Heitler process is a lot more efficient than the photo-pion process. Although relativistic 
shocks are likely capable of accelerating electrons to 7^ ~ 10'^ , Bethe-Heitler process could still be important if the number 
of electrons produced above the photosphere vastly outnumber the electrons expected to be in the GRB jet from simple 
charge neutrality. If the GRB has proton luminosity Lp given by Lp = ri~^L^, the comoving electron density is Ue = Up ~ 
2 X 10^r;~^I/-,,52rj^i?^g^ cm~"^. The number density of Bethe-Heitler produced electrons, ubh is 

R/r => = afarn^R/r. (51) 

Tie 

Since we want to restrict ourselves to above the photosphere, the optical depth GTTieR/T < 1 or 

^ <aj!^^lQ^riT2u;l(l + z)-\ (52) 

rie Ue 

where n-y is given by eq ((6)1. It is somewhat counter-intuitive, but the Bethe-Heitler process is likely to be most important in 
jets with lower baryon loading, i.e. when 77 is large. The Bethe-Heitler process could be important for 7^ <^ 10^ — especially 
if for some reason the Fermi mechanism is unable to accelerate electrons to this Lorentz factor in GRB shocks — but for these 
to account for the 100 MeV photons from GRBs via the synchrotron process requires a very large magnetic field and the 
luminosity carried by such a magnetic field would greatly exceed lO^'^ erg/s. Therefore, it seems that at best there might just 
be a small part of the parameter space for GRBs where the Bethe-Heitler mechanism could play an interesting role in the 
generation of prompt 7-ray radiation. 



4 PROTON SYNCHROTRON 

Massive particles have lower radiative losses than lighter particles, and therefore more massive particles are easier to accelerate 
in shocks. The maximum Lorentz factor that protons can attain is much larger than the maximum Lorentz factor of electrons. 
The maximum synchrotron photon energy from a shock-accelerated particle is given by requiring that the synchrotron energy 
radiated during one acceleration time (on the order of the Larmor time) is equal the energy gained in an acceleration cycle — 
half of the particle's energy, i.e. 'ymc/{qB) x AB^q'^j^ /(9m^c^) < 'ymc/2. The maximum photon energy for a source moving 
with a Lorentz factor F at redshift z is 



gpmc^ 50- 



167rq2(l + 2) l + z 



(m.) MeV. (53) 
\me/ 



While electron synchrotron radiation can only produce photons up to an energy ~ 50F/(1 + z) MeV, the proton synchrotron 
process can radiate photons up to 10^r/(l -I- 2:) GeV. For this reason, when photons of energies larger than what is allowed by 
elect ron synchrotron a re detected from a source, proton synchrotron is frequently suggested as a possible radiation mechanism 
Ce.g.lBottcher fc Dermer..l99a : Totani 1998; Aharonian 2000; Zh ang fc Meszaroa,2001. : Miicke et al.,.2003. : .Rcimcr ct al..i20oi : 



Razzague et al.ll2010l ) 



However, while the lower radiative efficiency of protons allows the protons to radiate at higher frequencies, it also means 
that the proton-synchrotron model requires more energy in the magnetic field to match an observed fiux. Because of this, we 
find that to match the typical observations of Fermi-LAT GRBs, either the energy requirements are prohibitive or the proton 
power-law distribution would have to begin at extremely high Lorentz factors. 

As before, we are considering protons with a power-law distribution dNp['yp) oc 7"'' d'yp if 7p ^ 7^. The proton injection 
frequency, Vi, is 



qBV% 



2-Kmpc{l + z) 



6.3 X 10~'"Br27ai + z)-' eV. (54) 



We define the cooling frequency, Vc, as the frequency where the synchrotron cooling time of the protons that radiate at Vc is 
equal to the dynamical time. The cooling time for protons is increased by a factor compared to the cooling time of 

electrons. The cooling frequency for proton synchrotron is 

z^c ~ 5 X 10^^B"^i?r5^ri(l + z)"^ eV. (55) 
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Efficiency 77 for Proton-Synchrotron Process 
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Figure 3. The maximum Tj = L^^/i^Lp -\- Lb) plotted in the Fl — F plane for various values of to match a flux of — 1 /^Jy at 100 
MeV. As in the previous figures, we fix 2 = 2, and dj^ = 4.9 X 10^* cm. This -q corresponds to the necessary luminosity in magnetic field 
and protons with energies greater than ^ 4r PeV to match a flux of = 1 fi]y at 100 MeV. The luminosity in protons with lower 
energies is likely comparable or much larger. Even when not considering the lower energy protons, the efficiency is very small for a vast 
majority of GRB parameter space. When r] is not shown either the place of emission is below the photosphere or it is opaque to radiation 
of 10 GeV due to 7 + 7 pair production. 



Since nearly all GRB observed in the Fermi-LAT band have a spectrum that can be fit by a single power law in the LAT 
band, we examine two possible spectral orderings: Vi ^ Vc ^ u and the slow cooling regime Vi ^ u ^ Vc- 

V must be above Vi to match the spectra of Fermi-LAT GRBs: v cannot produce GRB LAT emission because if the 
protons are cooled (uncooled), the spectrum is oc v'^ ^"^ {v^^^). These spectr a are harder than what is observed for most 



GRB, which have a typical high energy index < —1 ( Ackerman et al. 20121 '). Therefore, we take Vi v to agree with a 



typical GRB spectrum. The synchrotron flux at the peak of the spectrum [vi = min {vi, Uc)) is 

/. ~ 7B7V52r2(l + z)dl% fiJy, (56) 

P-i 

where N is the total number of protons radiating in a dynamical time. The flux scales as oc z/ 2 if ^ 1/ i/^ and as 

Below Vp, the flux of a typical GRB is constant, oc . Above Vp, the flux scales as fv oc v~^'^ . Since this break is larger 
than one half, it cannot be attributed to a cooling break. Therefore, in order to have oc below Vp, we require that both 
Ui and Uc lie above Up. Furthermore, the majority of LAT GRBs show a single power law above their peak, extending up to 
a maximum observed frequency, i/max, on the order of tens of GeV. We need to ensure that the proton synchrotron radiation 
does not add any spectral features in this energy range. Since we have already ruled out the fast-cooling regime, there are 
only two possibilities: the cooled case where, Ui ~ j/c ~ z^p ~ 1 MeV with a p ~ 2.4 and the uncooled case where Ui ~ Vp, 
Vc ^ I'max with p ~ 3.4. The cooled case can be ruled out because the energy required in the magnetic fleld is far too large. 
The uncooled case is considered in more detail in the following paragraphs. 

If we require that i/i ~ z^s (i-e. 100 MeV) and Vc ^ '-'max ~ 10 GeV, we can then place an upper bound on the magnetic fleld 
by requiring that Vc > i^max and a lower bound by requiring that protons will be able to be accelerated to high enough energies 
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to radiate at fmax- As a practical matter, these bounds do not affect our maximally efficient proton-synchrotron radiation 
calculation for the parameter range considered. We then minimize the total luminosity required in both the magnetic field 
and protons radiating at an observed frequency to match a typical observed flux of a few /ijy. 
The minimum Lorentz factor of the proton that radiates at is 

7, = 4 X 10«B-i/^r-i/^(l + z)'^^ul^'. (57) 

This Lorentz factor gives a proton luminosity Lp of 

Lp = 2r='7,7?-^mpc3iV ^ 5 x 10''\f,,^,jyB-''^^rl/^ R'r^^^^ {1 + 2)-'''^di,28 erg/s. (58) 

Given the magnetic field luminosity, Lb = 6 x 10'^'^ RIj^F?, erg/s, the total luminosity Lb + Lp will be minimized with 
respect to B when Lp = |Lfl, or when 

B ^ 10^ r-'^'RirfZy-s'"(^ + dZ, Gauss. (59) 

This magnetic field gives a proton luminosity 

T c w ir,52 j.4/7 -pl2/7 r,2/7 2/7/, , n-2/7,8/7 , ,„„^ 

Lp = -Lb ~ 5 X 10 f'jS^ Ri5 (1 + ^) ciL,28 erg/s. (60) 



The efficiency ri is plotted in figure [S] Note that Ls is not negligible as in figure [T] and [2l so in figure [S] = L^/{Lp + Lb)- 
The proton luminosity is much larger than the 7-ray luminosity for most of the allowed GRB parameter space. Additionally, 
proton synchrotron requires an unrealistically large 7;. Using eq (|59p . 

7, - 4 X io« i??rr-/^(i + .f'.rf-XdL%'- (61) 

It is unclear what physical process could produce a power law with such a high minimum Lorentz factor. The minimum Lorentz 
factor of a particle accelerated in relativistic shocks is approximately equal to the Lorentz factor of the shock front with respect 
to the unshocked fiuid, if every proton crossing the shock front is accelerated. The Lorentz factor can be proportionally larger if 
a small fraction of particles are accelerated and the remaining particles are "cold" downstream of the shock front. Considering 
that the Lorentz factor for GRB internal shocks is of order a few to perhaps a few tens, the typical proton Lorentz factor 
sho uld be ~ 10 — 10'^ (the large r value corresponds to when only 1 in 10^ protons are accelerated, as suggested by simulations, 
e.g. 



Sironi fc Spitkovskv ( 201ll ')). 7^ ~ 4 x 10^ is an unrealistically high injection Lorentz factor for relativistic shocks. If we 



set 7i to 10'', the proton synchrotron radiation would extend down to ~ 1 keV and would over produce in the GEM band. 
We can only decrease ji by a factor of 10 before over producing below the peak of the GRB spectrum. In summary, if proton 
synchrotron is to explain the observed LAT emission in GRBs, all of the protons must be accelerated to extremely high 
Lorentz factors very efficiently by some unknown mechanism. 

The expected neutrino fiux for the proton synchrotron model is estimated below. The total number of muon neutrinos is 
~ 2Tp^Np, where Tpj is the optical depth to photo-pion production, given in eq 0, and Np is calculated using eq ((56)) and 

imi): 

l\r o„in47r 1^-20/7 r)-l/7 rS/T -1/, , ^-13/7 -1/7,10/7 ,„„\ 

N^^2xl0 L^,52r2 R^^' /^;^jyi^p,6(l + 2) d^'^g. (62) 

Because the proton synchrotron radiation requires higher energy protons and larger magnetic fields compared to the photo-pion 
process, the pions produced will suffer larger radiative losses from synchrotron radiation before they decay. If the magnetic field 
is given by eq (|59p we find that the pions will be cooled significantly by synchrotron radiation when V^^^'^R^v] < 4/^'2jyi^8'"^(l+ 
z]^^"^ if'l\g. The neutrino fiux will peak at an observed energy of ~ ir7im^cV(l + z) ff the pions are uncooled. If the pions 
are cooled, the fiux will peak at the energy where pion cooling becomes important,or an energy of ~ ■jr7^.cooi?7i,rC^/(l + z) 
(7,,,ooi = \^^^B-^V2Rl^). 

{1.4 X lO'^ JZ]!3y^T i'^ + z)'^'''dl'^2s GeV if pions are not cooled 

(63) 
3.5 X l(fvf^''R:{'^f-%lv^^'\l + z)-^/''dl^l^GeV if pions are cooled 

Since the protons are not cooled by the synchrotron loss mechanism, the observed neutrino flux is calculated using eq H16|) 
and the dynamical time. The neutrino fiux, Fy, peaks at an energy given by eq (|63p and is 

{1.4 X 10"^ J^7,52r2 ^/^i?~//^fc'-^z/|''V*Cjy'^L%^(l + ^y'^^'^ GeV cm-^s"^ ff pions are not cooled 
(64) 
3.3 X 10"^ L^.52Tf^'' R'^'^-lv-^'"" fjjjydl]^l\l + 2)-"'"' GeV cm-^s"^ ff pions are cooled. 

As in Section [2TT7ll we estimate the neutrinos detected by IceCube per second of LAT emission from the proton synchrotron 
process: 
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( lO^"* L-y 52^2^^'' Ri^'^^^* f^^i^fyi^l^''^i^pld^'2j {^ + z)^'''^^'''"^ counts s^"*^ if the pions are not cooled 

^= ^/^^ 3^ ■ (65) 

(^5.6 X 10 * -^7,52r2^'^-Ri6^^^^^/^''Jjyi'8 ^^'^'^p.s'^L 28^(1 + •2^) counts s if the pions are cooled. 

Therefore, for a bright GRB detected by the Fermi- LAT with L^,j,2 ~ 10, /y,/iJy ~ 2, r2 ~ 9, 2: ~ 2 and fpx, ~ 1 and a 
duration of 10 seconds we find that the pions will be cooled if 7i < 2 x 10^'' cm. We expect ~ 4 x lO^^'-Rjg^^''^'' neutrinos of 
energy 2 x 10''R%'' GeV if i? > 2 x 10^-* cm and ~ 6 x lO-^E^g"''" neutrinos of energy 2.7 x 10''i?i^'' GeV if < 2 x 10"cm. 



5 SUMMARY AND DISCUSSION 

With a goal of understanding typical observed 100 MeV fluxes in bright FermirLAT GRBs during the prompt emission, we 

estimated the generation of photons by high-energy protons traveling through a shell of photons whose energy distribution 
is given by the Band function. We calculated the minimum energy in protons required to reproduce Fermi-LAT observations 
for the following hadronic processes: photo-pion, Bethe-Heitler pair production, and proton synchrotron. 

Unlike previous works, we specifically focused on the energy required for hadronic models to produce the >100 MeV 
photons seen in Fermi GRBs and how this requirement depends on GRB parameters. To provide additional physical insight 
into the energy requirements, we have provided both analytical estimates and more detailed numerical calculations. 

Wc find that photo-pion A-rcsonancc is much more efficient than Bethc-Hcitlcr pair production at producing high-energy 
electrons — so much so that Bethe-Heitler pair production can be ruled out as a mechanism for producing > 100 MeV photons 
observed by Fermi-LAT. The photo-pion process is capable of producing high energy electrons, but to match the FermirLAT 
flux at 100 MeV, the photo-pion process requires an energy in protons that is > 10'* times greater than isotropic energy in the 
7-ray photons. Since the Bethe-Heitler photo-pairs are produced more efficiently at low energies, Bethe-Heitler production 
will be the dominant process at low energies. These low-energy Bethe-Heitler electrons will have the same spectral index as 
the high-energy photo-pion electrons (assuming there isn't a cooling break). Therefore it is possible that both processes could 
add up to produce a single power-law deviation from the Band function that extends from low to high energies. This type of 
spectral feature has been observed in several Fermi GRBs. 

According to our calculations, proton synchrotron is capable of producing the Fermi-LAT GRB emission more efficiently 
than the other hadronic processes. The proton synchrotron could possible achieve efficiencies on the order of 1-10% for the 
brightest GRBs, if we assume that the minimum Lorentz factor for proton accelerated in shocks is extremely large ~ 2 x 10®. 
The minimum proton Lorentz factor of order 10^ is much larger than what is expected based on oxir current understanding 
of relativistic coUisionless shocks. Regardless of the mechanism accelerating the protons, these high energy protons with LF 
greater than 10® likely carry only a small fraction of the total energy carried by the protons. 

Wc also calculated the expected neutrino flux if the LAT emission is from the photo-pion process or proton synchrotron 
radiation. In the photo-pion process, for a bright LAT GRB with a Lorentz Factor of 900 and a duration of 10 seconds, we 
expect ~ .1 neutrinos detected by IceCube at an energy of ~ 10® GeV. Therefore, it may be possible to rule out photo- 
pion emission when the emission from multiple bursts is considered. For proton synchrotron radiation, the neutrino flux also 
depends on the emission radius, R. For a bright LAT GRB with a Lorentz factor of 900, an emission radius of 10^^ cm and 
duration of 10 seconds, we expect a ~ 4 x 10~^ neutrinos detected by IceCube at an energy ~ 2 x 10^ GeV. 

In summary, all the hadronic processes considered in this paper require signiflcantly more energy in protons than the 
observed energy in gamma-rays to reproduce the high-energy flux observed in Fermi-LAT GRBs. 
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